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ABSTRACT 

Recent observations have pointed out various anomalies in some multipoles 
(small £) of the cosmic microwave background (CMB). In this paper, it is proved 
that some of these anomalies could be explained in the framework of a modi- 
fied concordance model, in which, there is an appropriate distribution of vector 
perturbations with very large spatial scales. Vector modes are associated with di- 
vergenceless (vortical) velocity fields. Here, the generation of these modes is not 
studied in detail (it can be done "a posteriori"); on the contrary, we directly look 
for the distributions of these vector modes which lead to both alignments of the 
second and third multipoles and a planar octopole. A general three-dimensional 
(3D) superimposition of vector perturbations does not produce any alignment, 
but we have found rather general 2D superimpositions leading to anomalies sim- 
ilar to the observed ones; in these 2D cases, the angular velocity has the same 
direction at any point of an extended region and, moreover, this velocity has 
the same distribution in all the planes orthogonal to it. Differential rotations 
can be seen as particular cases, in which, the angular velocity only depends on 
the distance to a rotation axis. Our results strongly suggest that appropriate 
mixtures of scalar and vector modes with very large spatial scales could explain 
the observed CMB anomalies. 

Subject headings: cosmic microwave background — cosmology: theory — large-scale 
structure of universe 



1. Introduction 



The analysis of the data obtained by the Wilkinson Microwave Anisotropy Probe (WMAP) 
has pointed out some anomalies in the temperature distribution of the Cosmic Microwave 
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Background (CMB). These anomalies have not been explained in the framework of the con- 
cordance model, which is an inflationary flat universe with cold dark matter, dark energy, 
and reionization. For appropriate values of the involved parameters, this model explains 
most of the current cosmological observations, e.g., the magnitude-redshift relation satisfied 
by far supernovae, the statistical properties of galaxy surveys, and the CMB anisotropies; 
nevertheless, some aspects of these observations remain controversial. Among them, the 
WMAP anomalies deserve attention. Some of these anomalies could be due to unexpected 
systematic errors associated to foreground subtraction, galactic cuts, statistical analysis, and 
so on; however, other anomalies could be true effects requiring new physics. Future exper- 
iments as PLANCK should distinguish between physical effects and systematic errors. Let 
us now list the main anomalies: (i) the amplitude of the C 2 multipole is lower than it was 
expected, (ii) there is an asymmetry between the North and South ecliptic hemispheres, (iii) 
the multipole C 3 is too planar, and (iv) the multipoles C 2 and C 3 are too aligned. Other 
anomalies concerning £ > 3 multipoles have been also described. 

The importance of the anomaly (i) was initially overestimated. The probability assigned 
by ISpergel et al.l (120031 ) to the C 2 valu e obtained from the first year WMAP data was ~ 
1.5 x 10 -3 . Afterward, other authors (jEfstathioul 120031 ; iGaztanaga et al.l 120031 ; lEfstathiou 



2004 ; ISlosar et al.l I2004T) obtained grea ter probabilities by using different methods for data 



analysis. Finally, iHinshaw et al.l (120061 ) used the data from the first three years of the WMAP 
sky survey, plus appropriate statistical and foreground subtraction techniques, to conclude 
that the probability of the measured C 2 multipole is ~ 0.16. In conclusion, the observed value 
of C2 is currently considered small but compatible with the concordance model. Nevertheless, 
a lack of correlations at the largest angular scales appears to be s tatist ically significant in 



Spergel et al. 


(2003 


); 


Copi et al. 


(2007) 



The anomaly (ii) was studied in detail bv lEriksen et~aD (|2004al Jbh ; Iflansen et al~l (|2004al Jbh 



The hemispherical power asymmetry is nowadays considered substantial and rob ust, never- 



theless, more study is necessary to get definitive conclusions ( lEriksen et al.l 120071 ). 



Mathematical methods to quantify the alignment of C? and C3 as well as the planar 
charac ter of C 3 we r e depi cted by lde Oliveira-Costa et al.l (120041 ) [vectors n 2 , n 3 and parameter 
t] and I Copi et al.l (120041 ) [multipole vectors]. For the sake of simplicity, we have designed 
a code to compute n 2 , n 3 and t, whereas multipole vectors will be considered elsewhere. 
Vectors n 2 and n 3 make maximum quantity 



* = ^m 2 \a lm (n)\ 2 



[I) 



for £ = 2 and t = 3, respectively. In this last equation, quantities a^ m {n) are the spherical 
harmonic coefficients of the CMB map in a coordinate system where n coincides with the z- 
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axis. See Ide Oliveira-Costa et al.l (120041 ) for the e xplicit definition of parameter t. Anomalies 
(ni) and (iy) have been studied in many papers (jSchwarz et al.l 12004 ; iBielewicz et al.l 12004 ; 
Copi et al.l 120061 . 120071 ) . The planar shape of C3 has been confirmed in the bibliography, 
but this characteristic of the octopole is not very unlikely in the concordance model. More 
problematic is the strong alignment of C2 and C3. Some authors state that the multipole 
alignment is actually anomalous and also that the alignment extends up to I = 5. They 



suggest th e existence of a sy mmetry axis (|Land fe Magueiioll2005l . 120061 ; iBernui et al.l 120071 ; 
Choi 120071 ). Other authors (IRackic fc Schwarzl 120071 ) propose the existence of a preferred 
plane without rotational symmetry. This proposal suggests either a differential rotation 
viewed from an arbitrary point of the space, which should be outside the rotation axis, or 
a more complicated vortical motion with aligned angular velocities. Motions of this type 
-in extended regions- can be simulated with appropriate combinations of large scale vector 
modes. 

Finally, let us mention another CMB anomaly which has been found at sm aller angular 
scale s : a non Gaussian cold spot (~ 10° size) located in the South hemisphere (jVielva et al. 
2004 ICruz et al.lbood : iMartmez-Gonzalez et aubood ). 



An anisotropi c Bianchi Vlh, m o del has been recen tly considered (jJaffe et al.l 12005 



Bridges et al.ll2005l : Ijaffe et alil2006al lbl: [ Ghosh et al.ll2007l ) with the essential aim of explain- 
ing most of the above WMAP anomalies; however, the authors recognize that their model 
does not explain the observed acoustic peaks. Ot her authors have st udied the anisotropy 
produced by big voids with appropriate locations ( Inoue fc Silk 2006al lbl) to account for the 
mentioned anomalies. Motivated by the above considerations about symmetries and vorti- 
cal (divergenceless) motions, we propose here another possibility which may contribute to 
explain the large angular scale CMB structure: the existence of vector perturbations with 
large enough spatial scales. Here, the main features of the first £ multipoles produced by 
these vector modes are estimated in the framework of a concordance model. 

By using appropriate large scales, only their contribution to the first multipoles are 
significant and, consequently, there are no problems wi th the acoustic peaks. In the linear 
regime, scalar, vector and tensor modes (lBardeenlll980l ) do not couple among them; hence, 
vector modes can be separately studied. Vector modes are vortical peculiar velocity fields 
which do not appear in stan dard inflation; n evertheless, large scale vector modes may appear 
in brane - world cosm ologies (jMaartensll2000l ) and also in models with appropriate topological 
defects ( jBunnl 120021 ) . Whatever the origin of the vector modes may be, we are interested 
in their possible effects on the CMB when their amplitudes are appropriately normalized. 
Two effects pro d uced by the same type of large scale vector perturbations were studied in 
Morales fc Saezl (120071 ). Some basic aspects concerning these perturbations can be found in 
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this reference (hereafter, paper I). 

Our background is the so-called concordance cosmological model, with a reduced Hubble 
constant h = 10~ 2 H = 0.71 (where H is the Hubble constant in units of Km s~ x Mpc^ 1 ). 
The density parameters of vacuum energy and matter (baryonic plus dark) are = 0.73 
and Q m = 0.27, respectively. All these paramet ers are compatible w ith the analysis of the 



three first year WMAP data recently published (ISpergel et al.ll2006l ). 



Along this paper, Greek (Latin) indexes run from to 3 (1 to 3). Units are defined in 
such a way that c = k = 1 where c is the speed of light and k = 87rG/c 4 is the Einstein 
constant. The unit of length is the Megaparsec. Symbols a, r], and z stand for the scale 
factor, the conformal time, and the redshift, respectively. Whatever quantity A may be, Aq 
(A e ) stands for the value of A at present (CMB emission) time. Quantity is assumed to 
be unity. This choice is always possible in a flat background. 



2. CMB anisotropy 

The most general vector perturbation of a FRW universe is a fluctuation of the metric 
g a p, the four- velocity u a , and the traceless tensor E a p describing anisotropic stresses. 

In the absence of scalar and tensor perturbations, the gauge can be chosen in such a 
way that the line element reduces to 

ds 2 = a 2 (-d7] 2 + 2h i dx t dr] + 5ijdx l dx j ) , (2) 

where the perturbations of the g^ metric components have been written in the form hi = 
{hi,h 2 ,h 3 ) = h. 

From the matter four- velocity, u a = (u°,u), one defines the peculiar velocity v = u/vP. 

Finally, the condition E!y = is assumed along the paper, which means that there are 
no anisotropic stresses conditioning the evolution of v and h. 

Let us now calculate the CMB temperature contrast, AT/T, due to the above linear 
vector perturbations. From the equations of the null geodesies, the following formula can be 
easily obtained: 

AT - - - ■ ; P e 9hi , 

- v c0 -n-v ce -n-nn J ^T7 dr l ^ (3) 



where n is the unit vector in the observation direction and v c = v + h. In the case of 
linear vector modes, the integral can be calculated along radial null geodesies of the FRW 
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background, whose equations are fj = — r, 9 = cf) = (in terms of the spherical coordinates r, 
9 and (f) associated to x l ). The dots stand for derivatives with respect to the affine parameter. 

Functions h(r), f) and v(r), f) can be expande d in terms of an approp riate basis (the 
fundamental harmonic vectors, see iBardeenl (119801 ) and lHu fc White! ( 119971 )) to write 



— * 

h(f, rf) 



[B + (k, r/)e + (/t) + B~(k, rf)e~{Kj\ exp(ik ■ r) d 3 k 



(4) 



where k is the wavenumber vector, k is the unit vector kj k, and functions B + and B are 
the coefficients of the /i-expansion. A representation of vectors e + and e ~ is 



ef = {±kxk 3 /k-ik 2 )/aV2 , 

ef = (±k 2 k 3 /k + iki)/aV2 , 
ef = Tcr/kV2 , 



(5) 

(6) 
(7) 



where a — (kf + k^) 1 ^ 2 (see paper I). Hereafter, the following compact notation is used 

Vector v(rj, r) is expanded in the s ame way 



B+(k,r])e + (K) + B-(k,r])e-(K) = B 
using the coefficients v ± (k,r]). Quantities v' c 



— B^ are gauge invariant (IBardeen 



19801 ). Under the condition Eij = 0, quantities B ± (k,rj) decrease as a~ 2 in both the ra- 
diation dominated and the matter dominated eras (see paper I). Therefore, vector metric 
perturbations being significant at decoupling (the end of inflation) would be negligible to- 
day (at decoupling). During matter domination, the following formula holds: B ± {k,7]) = 
QHlVt m v^{k) / k 2 a 2 (rj) . Furthermore, functions are proportional to a -1 (constant) in the 
matter (radiation) dominated era. According to these comments, vector modes producing 
significant effects on the CMB should not freely evolve from the early universe. Either 
they are produced by exotic processes (brane-worlds, strings, and so on) close enough to 
recombination-decoupling or they must be maintained by some field producing an appro- 
priate ^ vector component (see paper I). Using the above expansions and evolution 
laws, the relative temperature variation due to the last term of Eq. ([3]) can be rewritten as 
follows: 



AT 
~T~ 



QH 2 Q m / -^-Mf) 
Jo a 2 {r) 



where F(r) 



F pq (r)n p n q and 



F„(tO = ~ifjp v±(k) e±W exp(z£ • f) cPk . 



(9) 



This last equation can be seen as a Fourier transform for each pair (p, q) of indexes. After 
these transforms are performed for appropriate boxes and resolutions, function F(r) and the 
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integral in Eq. (jSj) can be easily calculated for a set of directions defining a sk y CMB map. 



A HE ALPIx (Hierarchical Equal Area Isolatitude Pixelisation of the Sphere, see lGorski et al. 



(119991 ) ) pixelisation covering the sky with 3072 pixels is used in our simulations. 



Apart from the above CMB temperature effects, vector modes produce a rotation of the 
polarization direction (Skrotskii effect). As it was proved in paper I, the rotation angle is 

6iP = 3H 2 n m / [n-G(r}}, (10) 

Jo a (t) 



where 



For the line element ([21), the components of the angular velocity in momentum 



spac e 



are Wi = a 3 e i jkW : ' k . From this relation and the W mn components given by iBardeenl (119801 ). 
one easily gets -at first order- the following formulas: 

W 1 = iv±(4k z -efk 2 ) , (12) 

W 2 = tv c ± (e±k 1 -e±h) , (13) 

W 3 = iv±(e?k 2 -e±k 1 ) ; (14) 

hence, the equation hjWi = is identically satisfied. The resulting components only depend 
on the gauge invariant quantities = v^(k)/a and, consequently, the angular velocity is 
an appropriate vector field in order to discuss the properties of the vector modes and their 
superimpositions in a gauge invariant way (it is not the case of the peculiar velocity). 

Various appropriate choices of v^(k) are considered in next sections. In each case, the 
angular velocity and the resulting AT/T and 5ip maps are analyzed. For the AT/T maps, 
the angle formed by vectors H2 and 113 (giving the directions of the q uadrupole and octopole) 



and th e parameter t defining the planar character of the octopole (see lde Oliveira-Costa et al. 



( 120041 )) are calculated. 



3. CMB anisotropy produced by a single vector mode 

An unique vector mode k u is first considered. In this way, some ideas -which are basic in 
next sections to understand the CMB effects produced by superimpositions of these modes- 
are pointed out. For an unique mode, we can write: 

vtu(k) = vtJ(k - k u ) - (t£)**# + 4), (15) 
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where the complex numbers vf u = v cuR + iv^ uI fix the amplitude of the chosen mode and 
5{k — k u ) and 8(k + k u ) are Dirac-distributions. Equation ( Tl5l) implies the relation [f^(A;)]* = 
—v^ u (—k), which ensures that the components of the angular velocity in position space, as 
well as the temperature contrast AT/T and the Skrotskii rotation angle Sip are real numbers. 
Moreover, for an unique mode, the coordinate axis in momentum space can be chosen in 
such a way that k u = (k u i, 0, 0) with k u i = k u > and, then, Eqs. (jSJ-flTj) leads to: 



0; e 



i/V2; 



2 . 



(16) 



For the sake of simplicity in the notation, the x%, x 2 , and x% components of the angular 
velocity are hereafter denoted W x , W y , and W z , respectively The same notation is used for 
the components of any other vector in position space. From Eqs. ffT2l - (fT6]) one easily gets: 



W„ = 



w 9l 



k u V2 
k u V2 



\ V cuR V cuR, 



(V, 



cuR 



cuR. 



sin£ + (v 



cul 



V cul) COS Z 



J cul 



cul 



sin£ 



(17) 
(18) 
(19) 



where £ = k u -r '= k u r sin 6 cos and variables r, 9, and are spherical coordinates in position 
space. Analogously, From Eqs. (J9]), (|T5|) and ([TBI) one proves that the only non- vanishing 
components of F pq (r) are F 12 = W z /k^ and F 13 = —W y jk 2 u . As it follows from these relations 
and Eqs. (TT8l - (fl9]) , functions Fyi and F i3 depend on our choice of the complex numbers i>+ 
and v~ u . Once these numbers have been chosen, the integral of the r.h.s. of Eq. (jSJ) can be 
easily written as follows: 



AT 6V2~H%tt 



0"m 



n 



where A 



V cuR + V cuR, B 



k u 



cul 



J culi 



J cuR 



{An 2 + B n 3 )/ c + (C n 2 + D n 3 )I s 
C = -(vt uI + v: uI ),D 
i~ 2 {r) sin£ dr , 

i~ 2 (r) cos£ dr , 



J cuRi 



(20) 



(21) 



(22) 

n 1 = sin 9 cos 0, n 2 = sin 9 sin 0, and n 3 = cos#. The integrals (12T|) and (122]) are to be 
performed along each of the 3072 directions configuring our HEALPIx map from emission 
(r e ) to observation (r = 0). Afterward, the resulting map can be analyzed by using our 
numerical code specially designed to get n 2 , n 3 , and t. 

The value k u = 2tt/L u with L u = 4x 10 4 Mpc has been fixed and, then, for A = C = 
and B = D = 6.6 x 10~ 10 (mode [1]), vectors n 2 and n 3 appear to be perfectly aligned 
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in the direction (0,1,0) and the octopole is rather planar (t = 0.93). The total AT/T 
map is displayed in the top panel of Fig. [TJ The central and bottom panels of the same 
figure show the quadrupolar and octopolar components of this map. Figure [2] has the same 
structure but it corresponds to A — B — —C = D = 3.3 x 10~ 10 (mode [2]). In this last 
case, there is no alignment. The angle formed by the vectors ri<i = (0.037,0.706,0.707) and 
^3 = (—0.037, 0.706, —0.707) is very close to 90° and parameter t takes on the value t = 0.93 
(as in the first case). Other angles and t values appear for other choices of parameters 
A, B, C and D. These results strongly suggest that random superimpositions of arbitrary 
vector modes should not lead to aligned n2 and vectors. This fact is verified in next 
section by considering a rather general 3D superimposition. 

Finally, another type of vector modes (hereafter called w-modes) deserves particular 
attention (see § for applications). In this case, the coordinate axis in momentum space 
are chosen in such a way that k u = (k u i,k U 2,0) and, then, the conditions v^ u = v~ u = v cu 
are assumed. Complex number v cu can be put in the form v cu = \v cu \ cos/5 + i \v cu \ sin/3. 
Similarly, we can write k u \ = a u cos a and k U 2 = a u sin a. The effect of an unique w-mode is 
now considered. By performing the same kind of calculations as for previous isolated modes, 
one easily get: 

W x = W y = (23) 

v cu \ (cos (3 cos £ — sin (3 sin £) , (24) 

where £ = k u ■ r = a u r (n 1 cos a + n 2 sina). Furthermore, the associated temperature 
contrast is: 

AT V2y/2H2Q m \v. 



I c cos P — I s sin (3 



sin 2a 



2 2 -nV)— — + n 1 n 2 cos2a . (25) 



Thousands of maps, Mi, corresponding to different values of a and f3 have been obtained and 
analyzed. Parameters a u and \v cu \ have been fixed. Their values are a u = 7r/(2 x 10 4 ) Mpc' 1 
and \v cu \ = 3.3 x 10~ 10 . In Fig. [3] we display three of these maps corresponding to distinct 
w-modes; they are different, but the spots are always aligned along the equatorial zone and, 
consequently, as it has been verified, vectors n 2 and n^, are aligned along the direction (0, 0, 1) 
and, moreover, the octopole is very planar t ~ 0.94. This type of alignment and a high t 
value (planar octopole) appear in all the maps. Other values of a u and |t> c J have been 
considered with the same result. If we superimpose many of these maps, the vectors n 2 and 
H3 of the resulting map are not always aligned; in other words, any combination of linear 
modes lying in the plane (hi, k 2 ) with v^ = v~ = v c q does not lead to multipole alignments. 

This fact is not surprising taking into account that, for a given map, directions r?2 and 
n 3 maximize the quantity \& defined in Eq. ([TJ), which is nonlinear with respect to the 
ap m coefficients. Superimpositions of w-modes have been numerically analyzed in a simple 
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way, we have taken 1521 maps and, then, other 1521 maps Nj have been obtained 
according to the following formula: Nj = YH=i M- From the analysis of the Nj maps, the 
following conclusions have been obtained: (i) vectors n 2 and n 3 are aligned in the direction 
(0,0,1) for 1314 of these maps, which appear to have rather planar octopoles, (ii) in the 
remaining 207 cases, there are no alignments and the octopole is less planar. In Fig. HJ one 
of these cases is displayed, the spots of the bottom panel are not aligned in the equatorial 
zone {t = 0.46) and, then, the direction n 3 is not parallel to (0,0,1). Indeed, it has been 
numerically verified that these directions are almost orthogonal to (0,0,1) in most of the above 
207 cases. A theoretical proof of this orthogonality is not easy as a result of the particular 
form of the nonlinear definition of n 2 , n 3 , and t. In §[3J this type of vector modes (w-modes) 
will be superimposed to simulate differential rotations and other symmetric divergenceless 
motions. Then, the fraction of the superimpositions leading to n 2 and n% alignments will be 
experimentally found. 



4. 3D superimpositions of vector modes 

According to Eq. (jUJ), functions F pq (r) can be calculated by using the 3D Fast Fourier 
Transform (FFT). In order to do that, 512 3 cells are considered inside a big box with a size 
of 2 x 10 5 Mpc. In this way, the cell size is ~ 390 Mpc and, consequently, vector modes with 
spatial scales between 10 4 Mpc and 5 x 10 4 Mpc can be well described in the simulation. 
We can then calculate function F(r) to perform the integral in Eq. (jSj); in order to do that, 
the observer is placed at an arbitrary point located in the central part of the simulation box, 
where the Fourier transform is expected to be well calculated and, then, the integration is 
performed for each of the 3072 directions of the pixel centers. The variations of F(r) along 
the photon trajectories are smooth and, consequently, the integrations giving AT/T can be 
easily performed. Furthermore, in a central cube with 1.2 x 10 5 Mpc per edge (60 % of the 
box size in our simulations), we can place 5 3 observers uniformly distributed and separated 
by a distance of 3 x 10 4 Mpc. Then, quantity AT/T can be calculated for each of these 
observers; thus, from a given simulation, the information we obtain is greater than in the 
case of one unique observer located, e.g., at the box center. 

In this section, it is assumed (as in paper I) that vf 0R and vf 0I are four statistically 
independent Gaussian variables with vanishing mean, and also that each of these numbers 
has the same power spectrum. The form of this common spectrum is P(k) = Ak Uv , where 
n v is the spectral index of the vector modes and A is a normalization constant. Two values 
of the spectral index: n v — 1 and n v = 2 have been considered. The spatial scales are varied 
from 10 4 Mpc to 5 x 10 4 Mpc in all cases (only very small wavenumbers are considered). 
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Four realizations of this 3D random superimposition of vector modes have been performed 
for each spectrum and, then, 125 observers have been located as described above in each 
of the simulation boxes. Thus, 500 simulations of the CMB relative temperature variations 
obtained from the last term of Eq. ([3]) have been obtained. Moreover, the corresponding 
500 simulations of the term — v ce ■ n have been also found. In all cases, linearity conditions 
\h(r)\ <<< 1 and |u c (r)| <<< 1) have been verified using the relations: 

h( v , r) = -6Hfo m a~ 2 (ri) j ^1 e*(/c) exp(ik • r) d 3 k (26) 

and 

v^(r], r) = a' 1 ^) / vf (k) exp(ik ■ r) d 3 k. (27) 



The analysis of all these simulations have let to the following main results: (i) the term 
— v ce ■ n is negligible against the last term of Eq. Q. In Fig. we present one simulation of 
each of these terms for n v = 1. Numbers in the bottom panel (—v ce -n term) are much smaller 
than those of the top panel [last term of Eq. (J3J)]. Obviously, this comparison is independent 
on the spectrum normalization. We have verified that the average (C2) corresponding to the 
500 maps of the term — v ce ■ n is ~ 1/ 600 times smaller than the average calculated from 
Eqs. therefore, the term — v ce ■ n is hereafter neglected and our study is restricted to 

the maps obtained from Eqs. ©-Q; (ii) the angle 0^23 subtended by directions H2 and is 
smaller than 10° in nine of the 500 simulations for both spectral indexes: n v — 1 and n v = 2. 
These numbers are compatible w ith the 8.33 cases expec t ed fo r a random distribution of 



direction n 3 around a fixed n,2 (see lde Oliveira-Costa et al.l (120041 )): (iii) parameter t appears 



to be greater than 0.94 in 40 and 42 simulations in the cases n v — 1 and n v = 2, respectively. 
These numbers are to b e comp ared with 35, which is the corresponding number obtained by 



de Oliveira-Costa et al.l (120041 ) in the case of an isotropic Gaussian random field. All these 



considerations are independent on the normalization of the spectra. 

We can conclude that 3D random superimpositions of large scale vector models do not 
explain either the observed alignment of C2 and C3 («23 — 10°) or the unusually planar 
octopole (t ~ 0.94). However, the study of some 2D distributions of modes is worthwhile. 



5. 2D superimpositions of vector modes 

Special superimposition of vector modes are now considered. They are 2D superimposi- 
tions leading to divergenceless motions in long sized zones, which are hereafter called parallel 
vorticity regions (PVRs). In each of these regions there is a privileged direction. Inside the 
region, the angular velocity (describing the local vorticity there) is parallel to the privileged 
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direction everywhere. The a^-axis (hereafter z-axis) can be chosen to be parallel to the priv- 
ileged direction. Finally, the PVRs are assumed to be uniform along this axis in the sense 
that all the orthogonal planes are equivalent. In short, inside the PVRs, the components of 
the angular velocity are W x = 0, W y = and W z = W z {x x ,x 2 ). This configuration appears 
if functions vf {k) are chosen as follows: 

v+o(k)=v;o(k) = v c0 (k)5(6 k -^) . (28) 

In this equation, angle 9k is one of the spherical coordinates in momentum space (k and 
being the other two) and 5 stands for the Dirac distribution. By substituting the distributions 
in Eq. (|28|) into Eqs. (IT2l) - (fr^j) . the following relations are obtained in position space: 

W x — Wy — (29) 

W z (x 1} x 2 ) = V2 J v c0 (h, k 2 , 0) a^^^ dh dk 2 . (30) 
Analogously, from Eqs. (128j) and ([9]), the non- vanishing components of F pq appear to be: 

Fn = -F 22 = -V2 j ^ vaiki, k 2 , 0) e l ^ +k ^ dh dk 2 , (31) 

Fi 2 = V2 J M v^ki, k 2 , 0) e < k ^ +k ^ dk, dk 2 , (32) 

F 21 = -V2 J M Vc0 (h, k 2 , 0) e i{ - k ^ +k ^ dh dk 2 . (33) 
Finally, vector G(r) involved in Eqs. (|r0l) - (flT!) has the following components: 

G x = G y = (34) 

G z (x u x 2 ) = -V2 J v c0 (h, k 2 , 0) e i{klXl+h2X2) dh dk 2 . (35) 

As it follows from Eq. (1281) . our 2D superimpositions are combinations of the w-modes 
studied at the end of § [3] (k 3 = and v^ Q = v ~ Q = v c0 ) and, consequently, vectors n 2 and 
are expected to be either parallel or orthogonal (almost in all cases). The proportions 
between alignments and no alignments will be numerically obtained from the analysis of 
simulations. 



12 



5.1. 



Differential rotations 



A present angular velocity of the form W z = W z (p) is assumed, where p = \Jx\ + x\. 



This velocity describes a particular PVR, which could be interpreted as a big region under- 
going a differential rotation. The local vorticity only depends on the distance to the z-axis, 
which plays the role of the rotation axis. Then, from Eq. ( l30il one easily finds 



Function i> c o(&i> k 2 , 0) is calculated by using the last equation and, then, this function is 
substituted into Eqs. (|3T l) -( l33l) . to get the F pq components. It is also substituted into Eq. ( |35l) 
to obtain G z . All these functions only depend on the coordinates x\ and x 2 . They are 
easily extended inside a 3D cube (where photons move) taking into account that the planes 
orthogonal to the z-axis are indistinguishable. For example, in the case of function F 12 , 
its value at any point with coordinates (x 1 ,^ 2 ,^ 3 ) located inside the 3D cube would be 
Fi 2 (x 1 , x 2 , x 3 ) = Fi 2 (x 1 , x 2 , 0). These extended functions allow us to calculate either AT/T 
(from Eq. ([8])) or the polarization rotation angle 5ip (from Eq. (|T0l) ). These calculations can 
be performed for any observer located well inside the cube; in other words, for any observer 
whose last scattering surface is fully localized inside the cube. 

It is worthwhile to notice that, in the case of the rigid rotation of a big region, the 
angular velocity W z vanish in a certain gauge, in which the observer rotates with the region. 
In this gauge, Eq. (I36I) gives: v c (kx, k 2 , 0) = and, taking into account that this quantity 
is gauge invariant, it vanishes in any gauge; therefore, according to Eqs. (!3~T]) - (!35|) plus Eqs. 
([8]) and (jTUl) . quantities AT/T and Sip vanish. In short, there is no either CMB anisotropy 
or Skrotskii rotations associated to rigid rotations (the same is valid for rotations of the 
spatial coordinates in the absence of vector modes). These effects only appear in the case 
of differential rotations, which cannot be globally avoided by any rotation of the reference 
frame. 

Two functions W z = W z (p) have been used: the first one is 



where N\ is a normalization constant. The length m defines the spatial size of the PVR. 
The values m = 5 x 10 3 Mpc (case NI) and m = 3 x 10 3 Mpc (case Nil) have been tried. 
Evidently, the spatial scales involved in this differential rotation are very large. The second 
function is: 




(36) 




p < 2m 
p > 2m 



(37) 




P < P 
P > P 



max 



max 



(38) 



13 



quantities N2 and p max being the normalization constant and the parameter defining the 
spatial profile of the angular velocity, respectively. Two values of p max have been studied: 
Pmax = 6.8 x 10 3 Mpc (case CI) and p max = 4 x 10 3 Mpc (case CII). 

Once an angular velocity profile W z = W z (p) has been assumed (cases NI, Nil, CI, and 
CII), only two elements remain free: (i) the normalization constant, and (ii) the location 
of the observer in the simulation square. The square is that appropriate for the Fourier 
transforms in Eqs. (l3~Tj) - (l3"5]) . For the above profiles, a square size of 5 x 10 4 Mpc is used 
and, then, 81 observers are uniformly located in a central square of 2 x 10 4 Mpc size. The 
separation between neighboring observers is 2.5 x 10 3 Mpc; therefore, once parameter m 
(pmax) is fixed in the profile W z (W^), 81 simulations of AT/T and 5ip can be obtained as 
it has been described in the first paragraph of § 15.11 Each map corresponds to a localization 
of the observer characterized by its distance to the rotation axis (x\ = £2 = line). The 
analysis of the resulting HEALPIx maps has let to the following main conclusions: (1) the 
C2-C3 alignment is perfect for any of the above W z profiles and observers («23 = 0), (2) 
the inequality t > 0.94 also is satisfied in all cases. These results are encouraging. The 
proposed differential rotations plus appropriate large scale scalar modes could easily lead 
to the observed angle «23 — 10° and also to the parameter t ~ 0.94. Of course, the large 
scale vector modes under consideration should dominate against the scalar ones. Thus, the 
alignment produced by the differential rotation (vector modes) would not be hidden by the 
effects of standard scalar modes. The amplitude of the scalar perturbations contributing 
to small i multipoles (very large scales) should be smaller than those corresponding to the 
standard flat spectrum (compatible with the remaining observed Ci quantities). Either a 
certain cutoff or a damping of the scalar fluctuations would be necessary on very large scales. 
Details about the possible cutoff scale or the gradual damping are out of the scope of this 
paper; however, the general considerations of this paragraph are important to normalize the 
W z profiles. 

A few considerations about recent C MB observations are necessary before describing 



our normalization method. According to iHinshaw et al.l (120061 ) (WMAP three years data 
analysis), the CMB quadrupole is C^ MAP ~ 2.96 x 10" 11 whereas the octopole is C^ MAF ~ 
7.38 x 10~ n ; hence, if it is assumed that the contribution of scalar and vector modes to 
these multipoles are to be added (statistical independence of the scalar modes and the 
differential rotation) and, moreover, it is taken into account that the contribution of the 
vector modes must dominate (see previous paragraph), such a vector contribution should 
roughly satisfy the following conditions: (a) C2 must be a little smaller than 2.96 x 10~ n 
and, (b) 2C*2 < C3 < 3C2; hence, the following method is used to normalize in each of 
the cases NI, Nil, CI, and CII: in a first step, the C2 and C3 multipoles of the 81 maps 
are calculated for an arbitrary normalization and, then, the maps (observers) compatible 
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with condition (b) -which is independent on normalization- are found. The total number, 
Nb, of these maps is given -for each case- in Table [TJ Some of these maps correspond 
to observers located at the same distance from the rotation axis and, consequently, their 
normalizations are identical except for small numerical errors. This fact has been verified. 
The total number of distinct distances (observers), Nd, and the distances themselves, d or (i) 
with % : l,Nd, are also given in Table U) In a second step, the normalization constant is 
chosen to have C2 = 2.5 x 10~ n for each of the above observers and, then, the resulting 
octopoles, C$(i), are calculated and shown in Tabled] for % : 1, Nj. A number of different 
normalizations is thus obtained. Each of these normalizations is separately considered. The 
AT/T and Sip maps corresponding to one of the two observers of case NI (i — 1 in Table [1]) 
are displayed in Fig. [6j Top panel shows a AT/T map which seems to be clearly compatible 
with a planar octopole (estimated value: t ~ 0.9979) and a perfect alignment (a 2 3 = 0, with 
possible small errors due to the limited angular resolution of the HEALPIx maps). Bottom 
panel displays the corresponding 5ip map. Angles close to 0.1 degrees are reached in some 
directions, the angles are similar (a little smaller) than those obtained in paper I, which were 
estimated by using a rather arbitrary normalization. 
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Table 1. 2D simulations based on W z profiles. 



CASE 


N b a 




C 3 (l) x 10 11 


C 3 (2) x 10 11 


d or (l) X 10~ 3 


ci or (2) x 1(T 3 


x 10 9 


A wz (2) x 10 9 


NI 


12 


2 


7.23 


5.94 


7.9 Mpc 


9.0 Mpc 


0.99 


0.95 


Nil 


8 


1 


6.58 




7.9 Mpc 




2.51 




CI 


12 


2 


7.04 


5.33 


7.9 Mpc 


9.0 Mpc 


1.51 


1.41 


CII 


8 


1 


5.77 




7.9 Mpc 




3.94 





Note. — First column lists the four W z profiles defined in the text. In each case, 81 observers are uniformly distributed in the 
central part of the simulation box 

Q The number of observers whose CMB multipoles satisfy the relation 2C2 < C3 < 3C*2 is N b 

''Among the JVj, observers, there are N d ones which are actually different (they are located at distinct distances, d or , from the 
rotation axis) 

Note. — C*3(l) is the octopole (after normalization by the condition C2 = 2.5 X 10~ n ) of one of the Nd observers, whereas 
C3 (2) corresponds to the second of these observers (if it exists) . The same for d or and for the dimensionless ratio W z (p = 0) /Hq 
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After the above normalization method has been applied, any of the Nd normalizations 
corresponds to an observer (characterized by its distance to the rotation center) whose Ci 
and C 3 multipoles satisfy the following conditions: (i) they are appropriate to explain the 
values observed by the WMAP satellite with the help of a certain contribution due to scalar 
modes (to be estimated), (ii) these multipoles are fully aligned, and (hi) the octopole is very 
planar (t > 0.94). The distances from the observers to the rotation axis are different from 
zero (see Table [I]) and, consequently, these observers are not placed on the rotation axis but 
in another position, which is so much probable as any other position in the space. 

Normalizations lead to the values of the constants Ni and N2 involved in Eqs. fl37|) - fl38|) . 
from which, the dimensionless amplitude of the angular velocity profile A wz = W z (p = 0)/Hq 
can be found in each case. The resulting A wz values are given in TabledJfor the nor malizations 



inclu ded in it. They are a few times greater than the value 4.3 x 10 10 reported by (iJaffe et al. 



20051 ) in the framework of a fully different model. 



5.2. Statistical parallel vorticity fields 

In this section, a PVR region is simulated by using statistical methods. The compo- 
nents v c oR and v c qi of the complex numbers f c o(^i, 0) are generated as two statistically 
independent Gaussian variables with the same power spectrum and zero mean. The form of 
the spectrum is the same as in the 3D simulations; namely, -P(c) = Aa nv , and the chosen 
spectral indexes and spatial scales are also the same as in the 3D statistical realizations. 

Ten realizations of these 2D random superimposition of vector modes have been per- 
formed for each spectrum (n v = 1 and n v = 2) and, then, 81 observers have been uniformly 
located in the simulation square using the same method as in the 2D simulations with W z 
profiles; however, the sizes of the simulation square and the central square are 2 x 10 5 Mpc 
and 1.28 x 10 5 Mpc, respectively, and the distance between observers is 1.6 x 10 3 Mpc. 
Thus, 810 simulations of the CMB relative temperature variations produced by PVRs have 
been obtained. The corresponding 5ip maps have been also found. All these maps have been 
analyzed. Results from this analysis are now described; we begin with various conclusions 
which are independent on the spectrum normalizations: (a) the angle «23 is zero in the 
48.64% (48.4%) of the 810 simulations for n v = 1 (n v = 2), (/3) parameter t appears to be 
greater than 0.94 in the 18.64% (19.88%) of the simulations for n v = 1 {n v = 2), and (7) the 
conditions t > 0.94 and 0:23 = are simultaneously satisfied in the ~ 11% (~ 12%) of the 
simulations for n v — 1 (n v = 2). These last percentages can be found from Table [21 where 
the number of cases, n at , satisfying the two relations t > 0.94 and 023 = is given for each 
of the ten 2D realizations. We have counted these cases because, as it has been discussed in 
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§ [31 conditions t > 0.94 and 0^23 = do not seem to be independent and, consequently, the 
probability of the realizations satisfying the two relations is not a priori the product of the 
individual probabilities. 

The spectra are normalized as follows: first, all the simulations satisfying the conditions 
t > 0.94 and «23 = are normalized by the condition C2 = 2.5 x 10 -11 and, then, those 
of them satisfying the inequalities 2C*2 < C3 < 3C*2 are identified and counted. Their total 
number, n b s , is given in Table [2] for each of our ten 2D realization. It is worthwhile to 
notice that, each of the normalized simulations corresponds to one of the ten 2D statistical 
realizations and also to an observer located at a certain position in the simulation cube. 
Since there is no a rotation axis, coordinates x 1 and x 2 are both necessary to fix the observer 
position in the plane orthogonal to the vorticity direction of the PVR. 

For n v — 1 (n„ = 2), number n c & s appears to be zero in five (one) of our ten 2D 
statistical superimpositions of vector modes. In these five (one) cases, conditions t > 0.94 
and «23 = are satisfied (see Table [2]), but there are no observers measuring a quadrupole 
C 2 = 2.5 x 10 -11 and an octopole satisfying the relations 2C 2 < C 3 < 3C 2 . It is then easily 
calculated the probability of having at least an observer whose measurements satisfy the 
four conditions t > 0.94, 023 = 0, C2 = 2.5 x 10 -11 and, 2C*2 < C3 < 3C*2, namely, whose 
measurements may be compatible with current observations after introducing appropriate 
sub-dominant scalar modes. This probability is close to ~ 5.5% (~ 10.8%) for n v — 1 
(n v = 2). With these probabilities we cannot say that we live in a very special zone of the 
PVR, but in a reasonably probable one, which is equally probable than any other positions 
inside the PVR. 

For n v = 2 and the 2D realization number 9 of Tabled there are two observers (n b s = 2) 
whose measurements are compatible with the four above conditions. One of these observers, 
located at ~ 5 x 10 4 Mpc from the cube center, would measure t ~ 0.9631, a 2 3 = 0, 
C2 = 2.5 x 10 -11 , and C3 = 6.45 x 10~ n . The AT/T and Sip maps corresponding to this 
observer are shown in Fig. [71 Top panel displays the AT/T map, which looks like those 
compatible with a planar octopole and a perfect alignment. The corresponding Sip map is 
exhibited in the bottom panel. The largest angles - close to ~ 4.4 x 10~ 3 degrees- are much 
smaller (by a factor ~ 1/50) than those based on the normalization of paper I. Of course, 
these angles are too small to produce any currently significant 5-polarization of the CMB. 
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Table 2. Statistical 2D simulations. 



CASES (n v = 1) 


n a t a 




CASES (n v = 2) 


flat 


n obs 


1 


7 


i 


1 


11 


3 


2 


9 





2 


10 


1 


3 


6 


2 


3 


12 


2 


4 


9 


1 


4 


12 


2 


5 


9 





5 


11 


2 


6 


7 


3 


6 


10 


1 


7 


9 





7 


12 


2 


8 


4 





8 


12 


1 


9 


11 


1 


9 


12 


2 


10 


17 





10 


12 






Note. — Ten 2D statistical simulations corresponding to the spectral 
indexes n v = 1 and n v = 2 are numbered in the first and fourth columns, 
respectively 

"Number of observers which measure 023 = and t > 0.94 (among 81 
of them placed inside the simulation box) 

^Number of cases (among 81), in which measurements would be com- 
patible with conditions 023 = 0, t > 0.94, and 2C2 < C3 < 3C2 
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Finally, Fig. [8] shows a dimensionless quantity proportional to the present angular ve- 
locity W z . The represented zone is located inside the simulation square and centered in it. 
The normalization is the same as in Fig. [7] (same 2D simulation and observer). Red and blue 
spots correspond to regions which are rotating in opposite senses. A boundary with W z = 
separates them. The mean value of (W z /H)q is negligible by construction and the typical 
deviation is {\W z \ 2 ) l l 2 / Hq = 3. x 10~ 9 . Many realizations (as that of the Figure) have been 
considered to conclude that the typical value of (W z /H) is always a few times 10~ 9 . 



6. Discussion and conclusions 

Appropriate combinations of large scale vector perturbations have been introduced in 
the concordance model and, then, their effects on the CMB anisotropy have been studied in 
detail. Our main conclusions can be summarized as follows: 3D superimpositions of vector 
modes do not explain the CMB anomalies; however, some 2D superimpositions of these 
modes lead to good results. Two types of 2D simulations have been performed: one of them 
represents differential rotations of big regions and the other one leads to extended statistical 
PVRs. In these two cases there is a preferred direction of symmetry. It is the direction of 
the angular velocity, which is the same in any point of the perturbed region. In the first 
case, there is a symmetry around the rotation axis in the plane orthogonal to the preferred 
direction, however, statistical PVRs do not introduce such a rotational symmetry. 

Suitable differential rotations can explain the planar character of the octopole, its align- 
ment with the quadrupole, and the main part of the Ci and C3 values observed with WMAP. 
These facts are proved, in § 15. \\ for two different W z profiles. Polarization rotation angles 
8i[) close to 0.1 degrees are produced by these profiles. Other possible profiles could produce 
slightly greater angles. A sub-dominant contribution of large scale scalar modes could then 
account for a small part of the observed quadrupole and octopole, which would be comple- 
mentary of the part due to vector modes. These scalar modes could be also responsible for 
the observed angle 023 — 10°, which vanishes for pure differential rotations. The required 
scalar modes would destroy the rotational symmetry in the plane orthogonal to the axis of the 
differential rotation. Skrotskii rotations close to 0.1 degrees would produce a 5-polarization 
of the CMB, which could be marginally observable by future satellites (see paper I). 

For statistical PVRs, there is an appreciable probability of accounting for all the anoma- 
lies explained by differential rotations. This probability depends on the form of the assumed 
power spectrum and also on the interval of k values considered in the computations. The 
dependence on the spectral index has been pointed out by considering two distinct values 
n v — 1 and n v = 2 (see § 15.21) . The mentioned probability is greater in the case n v = 2 
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(~ 11%). Of course, a certain level of scalar modes is necessary (as in the case of differential 
rotations) in order to explain the observed angle 023 — 10°. Statistical PVRs lead to 8ip 
angles which are too small to produce significant levels of 5-polarization. Other intervals 
of spatial scales, and other power spectra could lead to higher probabilities for the explana- 
tion of anomalies and, perhaps, to greater Skrotskii rotations. In a certain k interval, the 
spectrum of vector modes could have any form (a power law is not required either by any 
theoretical prediction or by observational evidences). In a finite interval, e.g., between 10 4 
and 5 x 10 4 Mpc, the spectral index of a power spectrum is arbitrary, nevertheless, only some 
spectral indexes are admissible, as k tends to zero, to avoid divergences in some integrals 
(e.g., that of Eq. (13TD). 



In iRackic &: Schwarzl (120071 ). it is stated that, at high confidence, there is no any ro- 
tational symmetry of the CMB in the plane orthogonal to the symmetry axis. This fact is 
compatible with differential rotations by two reasons: (i) the mentioned rotational symmetry 
would be only observed from points placed on the rotation axis, whereas we are not located 
on this line with very high probability, and (ii) there may be either large scale sub-dominant 
scalar perturbations or deviations with respect to a perfect differential rotation and, obvi- 
ously, these perturbations and deviations could contribute to hide any rotational symmetry 
and also to explain the deviation from zero observed in the angle 023. 

The asymmetry of the the North and South ecliptic hemispheres is also compatible 
with our 2D superimpositions of vector modes. We predict two equivalent hemispheres, 
nevertheless, they are not separated by the ecliptic plane, but by the plane orthogonal to 
the angular velocity. Furthermore, in some slightly different scenarios, the equivalence of 
these two hemispheres could disappear. It occurs, e.g., if the last scattering surface of the 
observer is partially outside the PVR, which is particularly probable for PVRs which are not 
too extended in some direction. 

Solar system alignments would b e casual, as it seems natural in any cosmological expla- 



nation of the observed anomalies. See lChol (120071 ). All the theories of this type (including our 



proposal) would be ruled out by solutions of the CMB anomaly problem based on both, the 
ordinary spectrum of scalar perturbations, and a non cosmological component accounting 
for the observed statistical correlations with the local geometry of the solar system; how- 
ever, current observations and data analysis have not unveiled any component of this type 
accounting for the CMB anomalies. 

We have assumed very large spatial scales to alter only a few low-£ multipoles; never- 
theless, only vector modes have been considered. Why large scale scalar modes have not 
been tried? The main reasons are now pointed out. For the chosen spatial scales, combi- 
nations of modes should lead to very large almost-homogeneous regions. It occurs whatever 
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the nature of the perturbations may be. In the case of vector modes, the angular velocity 
will be almost- homogeneous in these regions and, consequently, it will have almost the same 
direction everywhere. These absolutely natural regions, which could have sizes comparable 
to that of the sphere bounded by the large scattering surface (for large enough spatial scales), 
are the PVRs we need to explain anomalies. In the case of scalar perturbations, the density 
contrast should be almost-constant in these large regions and, consequently, a cylindrical 
scalar inhomogeneity would be actually unlikely. Moreover, a flattened inhomogeneity does 
not seem likely as a result of the small scales required by the short thickness of the structure, 
which would be small enough to affect multipoles with too large £ values. Hence, symme- 
try axis and preferred planes seem to be rather improbable in the case of large scale scalar 
modes. Although these arguments are qualitative they strongly suggest the use of vector 
modes. 

Let us finish this paper with a list of a few open problems which should be addressed 
in the near future: (1) the origin and evolution laws of the vector modes (e.g. brane-worlds, 
strings, and so on) deserve particular attention. Only a consistent theory on these subjects 
could give answers to important questions as: in what a cosmological period (or periods) 
are generated the vector modes? How do they actually decay? How much probable are the 
PVRs? Are scalar and vector modes statistically independent? (2) Multipole components for 
i > 3 m ust be also analyz ed and compared with those extracted from WMAP data. Multipole 
vectors (ICopi et al. should be used in this extended study. (3) The proportions between 



large scale scalar and vector modes must be considered in more detail and, (4) deviations 
from the perfect parallelism assumed in our 2D superimpositions of vector modes could lead 
to interesting results (hemisphere asymmetry, «23 observed value, and so on). 

Large scale rotations are currently enigmatic (even for us), but the origin of the familiar 
cosmic expansion has kept unknown during a century. In both cases, rotations and expansion, 
rejection (acceptance) would be only justified by the disagreement (agreement) between 
predictions and observations (without prejudices). Although we have not a closed theory on 
the subject of this paper, results related with the CMB anomalies are actually encouraging 
and, consequently, more study is worthwhile. 
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Fig. 1. — Top panel shows the HEALPIx map of (AT/T) x 10 5 for the vector mode [1] (see 
text). Middle (bottom) panel displays the quadrupole (octopole) of this map. The alignment 
of C2 and C3 is evident. The octopole looks planar. Normalization is irrelevant. 




Fig. 2. — Same as in Fig. [T] for vector mode [2] (see text). There is no alignment in this case, 
but the octopole is visibly planar. Vectors r?2 and are almost orthogonal 
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Fig. 3. — Each panel shows the HEALPIx map of (AT/T) x 10 5 produced by a different 
w-mode (see text). The equatorial alignment and the planar character of the octopole are 
evident in the three panels. Normalization is irrelevant. 
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Fig. 4. — Same as in Fig CD for the superimposition of u>-modes Nn 2 e (see §EJ). Vector n 2 has 
the direction (0,0,1) and the octopole is not planar. Calculations indicate that vectors H2 and 
713 are orthogonal; hence, direction n 3 is contained in the equatorial plane. Normalization is 
irrelevant. 
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Fig. 5. — Both panels are HEALPIx maps of quantity (AT/T) x 10 5 . Top panel shows the 
part of this quantity due to the last term of Eq. ([3]), whereas the part obtained from the 
term — v ce ■ n is displayed in the bottom panel. The second part is much smaller than the 
first one. That is independent on the power spectrum normalization. 
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Fig. 6. — Top panel is a HEALPIx map of (AT/T) x 10 5 for case NI and observer i — 1 (see 
Tabled]). A strong C2-C3 alignment and a high t value are evident. The corresponding 5ip 
map is displayed in the bottom panel, where the angles are given in degrees. 
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Fig. 7. — Top: same as in Fig. [6] for one of the two observers measuring multipoles Ci and 
C*3 compatible with current observations in the realization number 9 of the case n v = 2 (see 
Tableland text). Bottom: the corresponding Skrotskii angles, Sip, are given in units of 
10~ 3 degrees. 




Fig. 8. — Map of the dimensionless quantity (W z /H) x 10 9 in the plane orthogonal to the 
angular velocity. The size of the represented square is 50 Mpc 



